Athermal Jamming vs Thermalized Glassiness in Sheared Prictionless Particles 
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Numerical simulations ol soft-core frictionless disks in two dimensions are carried out to study 
behavior of a simple liquid as a function of thermal temperature T, packing fraction <f>, and uniform 
applied shear strain rate 7. Inferring the hard-core limit from our soft-core results, we find that it 
depends on the two parameters <j> and T/7. T/7 — ¥ defines the athermal limit in which a shear 
driven jamming transition occurs at a well defined <f>j. T/7 — > 00 defines the thermalized limit 
where an equilibrium glass transition may take place at a (pa. This conclusion argues that athermal 
jamming and equilibrium glassy behavior are not controlled by the same critical point. Preliminary 
results suggest (f>G < 4>j- 

PACS numbers: 64.70.Q-, 45.70.-n, 64.60.-i 
Keywords: 



Seemingly disparate physical systems are found to un- 
dergo a transition from a flowing liquid-like state to a 
rigid but disordered solid state upon varying some ex- 
ternal control parameter pQ. Athermal granular parti- 
cles undergo a sharp jamming transition as the packing 
fraction <f> is increased above a critical value (f>j. Simple 
liquids may freeze into a disordered glass as temperature 
T is decreased below the glass transition T g . Foams may 
cease flowing and show elastic response as the applied 
shear stress a is decreased below the yield stress ay . An 
interesting conjecture by Liu and Nagel and co-workers 

[3] attempted to unite these phenomena into a single 
jamming phase diagram, with the three axes <p, T, and 
a. 

In the equilibrium qb — T plane, it may be expected 
that as 4> decreases, the glass transition T g {<f>) decreases 
and vanishes at a particular value <pc- In their original 
jamming phase diagram, Liu and Nagel proposed that 
this 4>g was the same value as the </>j that locates the 
athermal jamming transition. The point T = 0, <f> = <p,j, 
a = thus locates a special critical point J that con- 
trols not only athermal jamming, but also finite temper- 
ature glassy behavior. If correct, such a scenario could 
offer a new approach towards understanding equilibrium 
glassy behavior, by studying the effect of thermal fluctu- 
ations about the presumably simpler and better under- 
stood jamming point J. 

However, for protocols restricted to the <f> — T plane, 
it is now generally accepted from theoretical mean-field 
calculations [3HZ] and demonstrated by numerical simu- 
lations [8HIU] that the location of the T = jamming 
transition is in general protocol dependent. A contin- 
uous range of values for (j>j may be found, depending 
upon the particular ensemble of initial states, as well as 
on the rates of compression or cooling. Rapid compres- 
sion or cooling from random configurations leads to a 
lower bound, often associated with random close pack- 
ing. It remains unclear if a well defined nontrivial upper 
bound for (pj exists for vanishingly slow compression or 
cooling, or if crystallization must in principle occur for 



sufficiently slow rates [TT]. In contrast, a true equilib- 
rium glass transition T g should be protocol independent. 
The connection between athermal jamming and an equi- 
librium glass transition is therefore subtle, if one stays in 
the 4> — T plane. 

Here we address this problem by looking along the 
third axis of the Liu-Nagel phase diagram, considering 
the behavior of systems undergoing uniform steady state 
shear, with a fixed applied shear strain rate 7. Consid- 
ering a simple model of overdamped, frictionless disks 
in two dimensions (2D), we have earlier shown [9] that 
athermal steady state shearing defines a statistical en- 
semble of states, such that there is a uniquely defined 
shear-driven jamming transition cj>j in the limit of van- 
ishingly small strain rate, 7 — > 0, independent of the 
initial configuration. We now extend these investigations 
to finite T. We argue that the distinction between shear 
driven athermal jamming and thermalized glassy behav- 
ior may be considered within a strictly hard-core model. 
Using simple dimensional arguments, and by inferring the 
hard-core limit from simulations of soft-core particles, we 
show that athermal jamming and thermal glassy behavior 
are described by opposite limits of a key control param- 
eter of the hard-core system. We thus conclude there is 
no reason to expect that these two phenomena are con- 
trolled by the same physical processes, nor any reason to 
expect 4>g = 4>J- 

Our model is one that has been extensively studied at 
T = O0. We use a bidisperse mixture of frictionless soft- 
core disks in 2D, with equal numbers of large and small 
particles with diameter ratio di/d s = 1.4. The soft-core 
interaction between overlapping particles i and j is har- 
monic, V(rij) = eVfaj) = ^eSfj, where <% = (l-r^/dy) 
is the relative particle overlap, is the particles center 
to center distance and dij is the sum of their radii. We 
use Durian's bubble dynamics [12] of overdamped parti- 
cles with a viscous dissipation with respect to an imposed 
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average linear shear velocity flow in the x direction, 
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where A is the viscous damping coefficient. Temperature 
is modeled by a random Langevin thermal force Ci satis- 
fying, 



(C) = o , (c i (t)C j (*')) = 2Ar^(t-Oi 



(2) 



with I the identity tensor. We use Lees-Edwards bound- 
ary conditions |13| to apply a uniform shear strain rate 7 
to a box of length and height L containing N particles. 

Defining the time constant to = Xd 2 /T, we can cast 
the above equation of motion into dimensionless form by 
dividing each term by T/d s to get, 
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with dimensionless variables fj = ri/d s , t = t/t , and 
noise Ci satisfying the correlation (Ci{t)Cj(i')) = 2SijS(t— 
r)I. In the hard-core limit, e/T — > 00, the first term on 
the right hand side acts to provide an excluded volume 
effect, preventing particle overlaps but introducing no en- 
ergy or time scale. Behavior in this sheared hard-core 
limit is thus entirely determined by two dimensionless 
parameters: the packing fraction 4>, and the Peclet num- 
ber 7^0 oc 7/T (for later comparisons, we will find it con- 
venient to phrase our discussion in terms of the inverse 
Pcclct number oc T/7). This immediately yields one of 
our main conclusions. The effects of temperature and 
shear on the hard-core system enter only via the combi- 
nation T/7. The limit in which athermal jamming takes 
place corresponds to T —> first, followed by 7 — s- 0, 
i.e. the limit T/7 — > 0. The limit in which thermalized 
glassy behavior takes place corresponds to 7 — > first, 
followed by T — > 0, i.e. the limit T/7 — > 00. As athermal 
jamming and thermalized glassy behavior occur at the 
extreme opposite limits of the control parameter T/7, 
there is no reason to expect that jamming and the glass 
transition share a common physical mechanism, or that 
they occur at the same value of the packing fraction 4>. 

We now consider the shear viscosity r\ = er/7. In terms 
of the dimensionless variables defined above, the dimen- 
sionless viscosity 77, with a well-defined hard-core limit, 
is 77 — ri/Xd 2 s . By measuring the viscosity r/(cj),'y 1 T) 
of soft-core particles, we will infer the hard-core limit 
f/hc(0;T/7). The athermal jamming transition is then 
defined by the cj)j where ?jhc(4>J, 0) — > 00, while the 
hard-core glass transition is defined by the 4>g where 

%c(</>G,°o) -> oo- 

Our simulations are at a fixed packing fraction cj>, using 
N = 65536 particles so that finite size effects are negligi- 
ble for the parameter ranges utilized in this work. The 
elastic part of the stress tensor p a p is computed from the 
contact forces in the usual way [3J; the elastic part of 



the pressure is p = (p xx +p yy )/2 and the shear stress 
is a = p xy . Henceforth we measure length in units such 
that d s — 1, energy in units such that e = 1, and time in 
units such that Xd 2 /e — 1. In these units we have fj = r). 
Similar shear driven simulations at finite T have been 
carried out for underdamped particles by others |14H16j . 

We consider first a moderately dense value oi<j> = 0.72, 
but still well below the jamming <fij w 0.843 [T7]. In 
Fig. [IJ, we plot our results for 77 vs T for several differ- 
ent values of 7. We also show the linear response ?7gk, 
computed using the Green-Kubo formulaJTB] applied to 
equilibrium (7 = 0) simulations of Eq. ([T]) [T5]. We see 
that at high T all the data collapse to a common curve; 
•q is independent of 7 indicating a linear rheology where 
thermal fluctuations dominate over shear induced fluctu- 
ations. As T decreases, 77gk increases and appears to sat- 
urate to a fixed value, which is just the hard-core value of 
viscosity in thermal equilibrium, rfrceq- As T decreases at 
finite 7, however, 77 increases to a peak value 77poak(</ ) , 7) 
at a temperature T poa k(0, 7), and then decreases to a fi- 
nite value as T — > 0. As 7 decreases, 77 pca k saturates to 
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FIG. 1: Shear viscosity 77 vs (a) T and (b) T/7, at fixed 
packing fraction (j> — 0.72 below jamming, for several different 
shear strain rates 7. In (a) solid dots are the linear response 
»7GK computed from the Green-Kubo formula in equilibrium. 
In (b) the horizontal dashed line is the hard-core equilibrium 
limit 77hco q ; the horizontal dotted line is the athermal value r/o 
at T = 0. 



On the low T side of 77 pca k, the rheology is highly non- 
linear (77 varying with 7) until converging to a common 
value as T — > 0. This T = value is just the athermal 
viscosity 770. This low T side of 77 pea k can be better un- 
derstood by plotting 77 vs T/7, as shown in Fig. [T]d. We 
now see that the data below 77 pea k for different 7 collapse 
to a common curve. This common curve, as 7 — ¥ 0, rep- 
resents the hard-core limit predicted by Eq. ([3]), which 
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we see is a smooth sigmoidal shaped curve increasing 
monotonically from 770 in the athermal T/7 — > limit, 
to ?7hccq hi the thermalized T/7 — > 00 limit. Data at 
finite 7 to the right of ?7poak(7), that fall below this lim- 
iting 7 — > curve, represent the parameter region where 
soft-core effects are important. From Fig. [T^, we con- 
clude that as 7 — > 0, the hard-core region of the system 
gets pushed down to T — > 0. The thermalized limit, 
lim-r^o [lim-y_>.o v] = ^hccq, becomes singularly decoupled 
from the athermal limit, lim^o [liniT^o v] = ? /o- 

The non-monotonic behavior of t](T) at fixed 7 can be 
understood as due to competing effects of thermal fluc- 
tuations on the hard-core vs the soft-core regions of the 
system behavior. As T increases in the hard-core region, 
the forces associated with collisions increase and hence 
pressure p increases. Since particles cannot pass through 
each other, it is difficult for shear stress to relax and so 
as p increases, so does a and hence 77. As T increases 
further, however, one enters the soft-core region where 
there is enough thermal energy for particles to press into 
each other. Particles may now squeeze past each other, 
allowing for more rapid relaxation of shear stress, with a 
decrease in a and hence 77. 

It is interesting to consider how other quantities de- 
pend on the hard-core parameter T/7. In Fig. [2^, we 
plot the inverse reduced pressure nT/p. with n = N/L 2 
the density of particles, vs T/7. In the thermalized limit 
T/7 — > 00 we see that the curves, as 7 — > 0, approach the 
hard-core equilibrium value, as we have computed inde- 
pendently from pair correlations |19j evaluated in equilib- 
rium hard-core Monte Carlo simulations. In the athermal 
limit T/7 — > 0, nT/p — >• as is expected. In Fig. [2]d we 
plot the stress anisotropy a jp vs T/7. In the thermalized 
T/7 — > 00 limit we find a jp — > as expected; as 7 — > 0, 
the thermalized system maintains a finite pressure but 
no shear stress. In the athermal T/7 — > limit, a/p ap- 
proaches a finite value, as has been observed earlier [5D] ; 
in the athermal limit, both a and p vanish as 7 — > 0, but 
do so in a way that their ratio becomes constant. This 
appearance of anisotropy as one moves from the thermal- 
ize to the athermal limit has recently been observed in 
experiments on colloidal particles |21) . 
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FIG. 2: a) Inverse reduced pressure nT/p and b) stress 
anisotropy a/p vs T/7, at packing fraction <f> — 0.72 for sev- 
eral different shear strain rates 7. Dashed horizontal line in 
(a) is the hard-core equilibrium value of nT/p. 



Returning to viscosity 77, we now consider behavior at 
higher packing fractions 4>. In Fig. [3] we show results for 
77 for different values of 7 at the higher value of </> = 0.80. 
As in Fig. [T] we see that the data collapses on the high 
T side of 77 pea k when plotted vs T, but collapses on the 
low T side of rj pca ^ when plotted vs T/7. Unlike Fig. [l] 
however, we see that 77 pea k continues to increase as 7 
decreases, with no sign yet of saturating. Our data in 
this high T/7 limit is not at sufficiently small 7 to have 
reached the hard-core limit. Equilibrium simulations at 
this high tfi cannot be sufficiently equilibrated to directly 
compute ?7hcoq- 




FIG. 3: Shear viscosity rj vs (a) T and (b) T/7, at fixed 
packing fraction cj> = 0.80 closer to the jamming 4>j = 0.843, 
and several different shear strain rates 7. 



That the hard-core viscosity 77h c (</ ) , T/7) appears to 
be monotonically increasing as T/7 increases, and that, 
comparing Figs. [T] and [3j ryhcoq appears to be increas- 
ing much more rapidly than 770 as (f> increases, suggests 
a <f>G < 4>j = 0.843, in agreement with recent equilib- 
rium simulations of soft-core particles in three dimensions 
[22| 123] . In Fig. [4 1 we plot 77 pC ak vs 7 for several differ- 
ent values of <f>, while in Fig 4Jd we plot T pea k- Since, as 
7^0, 77 poak 77hccq, ?7poak should stay finite for <f> < 4> G , 
while 77 pca k — > 00 for <j) > <j)Q. Looking at the raw data 
in Fig. [4^,, such a change in behavior appears to happen 
at 4>* 0.80. T pea k shows a similar marked change in 
behavior at the same <fr* , with T pea k — ¥ as 7 — > for 
</><</>*, while T p e a k is decreasing much more slowly, and 
perhaps saturating to a finite value, for <fi > <p* . 

We can try to quantify this behavior with a critical 
scaling analysis. Assuming the usual algebraic scaling of 
a continuous critical point at <pQ 1 we expect 77 to satisfy 
a scaling relation, 
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FIG. 4: a) Peak value of viscosity 77 pea k, and (b) location of 
the peak viscosity T pGa k, vs strain rate 7 for different packing 
fractions 6. 



FIG. 5: a) (a) Scaling collapse of ?7 pca k according to the scaling 
relation of Eq. |6j. (b) Mean squared displacement (Ar 2 ) vs 
MC passes for different <j>. Inset shows diffusion constant D 
vs <j>. 



where 5<f> = 4>— 4>g and b is an arbitrary length rescaling 
factor. The peak in 77 as T varies for fixed (f> and 7 will 
occur when the scaling function /(x, y, z) has a maximum 
at some z vea ^(x,y). This leads to the scaling equation, 

V P cMi) = bP/»g(6^,jb z ) (5) 

with g(x,y) = f(x,y,z peak (x,y)). Choosing b = j' 1 ^ 2 
then gives 

r] pea ^ /zu = h{5^/j l 'n (6) 

with h{x) = g(x, 1). Expanding h(x) as a polynomial 
in x, we fit our data in Fig. [4ji to the scaling form of 
Eq. keeping (3, zv, (f>o, and the polynomial coef- 
ficients as free fitting parameters. Restricting to data 
with 7 < 10~ 6 our fit yields the data collapse shown 
in Fig. [5}i, with fitted values G « 0.796, (3 w 2.7, 
zv w 5.1. Choosing 6 = |<5</>|~ 1 ' in Eq. Q to get 
%cak = |5</>| _/3 g(±l, 7/|J0| ZI/ ), we see that j3 is just the 
exponent that describes the algebraic divergence of the 
thermalized hard-core viscosity as <fi — > <\>q. 

Although our fit in Fig. [5ji appears to look good, we 
find that the values of our fitted parameters arc some- 
what sensitive to the range of data that we use. As we 
restrict data to a smaller window of 7 < 7 max the fitted 
value of 4>g appears to increase slightly, while the x 2 /dof 
decreases. For our smallest window with j max = 10~ 6 
we find x 2 /dof ~ 6, suggesting a fit that is still less than 
ideal. Such behavior could be due to the effect of cor- 
rections to scaling, as we have previously shown to be 
important at the athermal jamming transition [17) . 

As an alternative approach to locating the thermalized 
4>Gt we have carried out independent equilibrium (7 = 0) 
Monte Carlo (MC) simulations of hard-core disks. At 
each step of the simulation a particle is picked at ran- 
dom and displaced a distance Sr with Sx, Sy chosen from 
a uniform distribution on the interval (—0.05, 0.05). If no 
particle overlap results, the move is accepted, otherwise 
it is rejected. N such steps constitutes one MC pass and 



represents one unit of time. With this MC dynamics we 
compute the single particle mean square displacement, 
(Ar 2 (t)) (1/N) Eili(|ri(t) - r,(0)| 2 ). In Fig. §3 we 
show our results for several different 0. We see that (Ar 2 ) 
continues increasing with time for <fi < 0.805, suggest- 
ing a finite diffusion constant D, however D is rapidly 
decreasing as <j> increases (see inset to Fig. [5}d). These 
results thus suggest a <j>G ^ 0.805. 

The discrepancy in the estimate of 4>g from our viscos- 
ity vs our diffusion data could be due to several factors: 
(i) the neglect of corrections to scaling in our analysis of 
(ii) it may be that the scaling shown in Fig. [HJi re- 
flects a mode coupling transition [Ml |2"5"] that gets cut off 
by some other physical mechanism on longer time scales 
as 7 decreases; (iii) diffusion at the higher (j> seems corre- 
lated with increased particle segregation [T7] , and so may 
be reflecting an approach to a true phase separated equi- 
librium [TT] rather than a metastable glassy state. We 
note that recent shearing simulations in three dimensions 
[26] similarly suggest a lower value for <pQ from viscosity 
measurements than was previously found from relaxation 
time measurements of both hard-core and soft-core 
[22j particles in equilibrium. The precise value of <pa an d 
a complete understanding of the thermal glass transition 
thus remains for future work. 

To conclude, we have demonstrated that in the hard- 
core limit, a system of overdamped steady-state sheared 
particles is characterized not only by the packing fraction 
(/>, but by an additional control parameter T/7. The limit 
T/7 — > corresponds to the athermal limit in which a 
sharp jamming transition at <f>j is clearly observed. The 
limit T/7 — > co corresponds to the thermalized limit 
where equilibrium glassy behavior is observed. That the 
athermal and the thermalized regions are at opposite lim- 
iting values of T/7 argues that there is no reason to ex- 
pect athermal jamming and thermalized glassy behavior 
to be controlled by the same critical point. Our results 
find behavior consistent with a 4>g < <pj, however we can 
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not yet be conclusive about the exact value of <j>a- As 
we were finishing this work, we learned of recent work by 
Ikeda et al. [55] who have reached similar conclusions. 
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